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♦ Introduction 

• Motivation 

• Near-field and far-field plume flow fields 

• Scaling Theory 

♦ Phoenix Mars spacecraft 

• Experimental 

• Numerical 

• Data comparison/Flow Physics 

♦ Mars Science Lab (MSL) Descent Stage spacecraft 

• Experimental 

• Numerical 

• Data comparison/Flow Physics 


♦ Conclusion 


♦ Plume-surface interactions due to spacecraft iandings 

• Spacecraft stability and survival 

- Moments/Torques 

- Updraft plumes 

- Plume induced heating 

• Cratering & dust lifting 

• Implications for manned and large payload landings to Mars, 

moon, asteroids and other planetary bodies 

♦ Theory, test data and numerical simulations were used to 
characterize the complex plume impingement physics 
and identify the environments observed due to 
spacecraft landings 



JPL 


1998 Mars Polar Landar Failure Report 
(WhetSBl el al, 2000) 

- MPL Project failed to conduct studies on plume-surface interactions 

- Recommended these investigations for future powered descent 
landing missions to Mars 




These NASA spaceflight projects deemed it 
necessary to conduct these investigations thru 
University research partnership. 

Phoenix: Pulsed-modulated descent system 
(Rocket Engine Module - REM) 


2007 


2011 


Mars Science Lab: Sky-crane throttled landing 
system (Mars Landing Engine - MLE) 


Last detailed study was completed in 1973 for Viking 1 
and 2. 
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Underexpanded Supersonic Jet 



Planar Laser Induced Fluorescence Imaging 


Overexpanded Supersonic Jet 



Inmann et a!., 2009 




Far-field flow/ 

Impingement zone - TEST DATA 

Important flow structures with implications 
to cratering, acoustics and spacecraft 
dynamics during descent 



Stagnation eubbia Tail Shock 



Lamont and Hunt, 1976 
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Normalized parameters 


Schlichting and Gersten, 2001 
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e = jet expansion ratio 


Free surface boundary condition 
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Nondimensionai numbers that satisfy dynamic simiiarity 
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Phoenix REM 




Vi scaie 

fuii-scaie 

Vs scaie 

fuii-scaie 

Hypersonic similarity 

k 

14.8 

14.0 

12.7 

11.4 

Jet expansion ratio 

e 

2.9 -2.1 - exp 

6.8 -2.2 -fit 

~4.4 - exp 

3.8 - fit 



3.5 - num 

6.8 - 4.1 - num 

4.5 - num 

4.7 - num 

Reynolds Number 

Re 

24.5-14.7x105 _exp 

8.4-5.0X 105-flt| 

12.7X 105 

3.4 X 105 



23.3 X 105 - num 




Mach Number 

Me 

5.14 

5.08 

4.77 

4.67 

Strouhal Number 

St 

0 

0 

4.4 X 10^ 

3.3 X 10^ 


Table 1. Scaling parameters; exp-experiment; num-numerical simulation; fIt-spaceflight conditions 




Courtesy of NASA/JPL/Lockheed Martin 
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Experimental Methodology 


Absolute Pressure Sensor (1) 



-V 2 scale Phoenix nozzle (MR-107) 

-N 2 test gas 
-10 Hz pulsing 
-Mars atmosphere 

University of Michigan Thermai Vacuum 
Chamber 


Numerical Methodology 


-Two Navier-Stokes computational solvers were 
used for modeling and analyses 

- ANSYS FLUENT 

-3-D & axisymmetric density based solver 
-Transient RANS 

-Time step - 1 us - explicit marching 
-Turbulent (RNG) model 
-Adaptive meshing for resolving shocks 
-Grid independence 

-2*^^ order upwind discretization scheme 
-2 million unstructured grid cells 

- Aerosoft GASP 

-3-D density based solver 

-Transient & steady-state RANS eqns solved 

-Van Leer flux splitting 

-Laminar 

-Dual implicit time stepping 
-Single species - frozen flow 
-Grid independence 
-4 million unstructured grid cells 

Plemmons, Mehta et al, 2008 



P / Pc-max 


MARS ATMOSPHERE ~ 700 Pa 


Temporal Profile Spatial Profile 

Observed large transient overpressures during Observe a monotonic rise in ground pressure 

engine start-up and shut down followed by a drop in centerline pressure 


Overpressure 



Back pressure 


Plemmons, Mehta et al, 2008 
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meters 


Mechanism deduced from experimental 
measurements, transient numerical simulations 
and theory. 


MARS ATMOSPHERE ~ 700 Pa 


Plate (recovery) shock formation 


Pc -1200 kPa 
0 deg cant 
h/de = 8.4 
N 2 test gas 

Axisymmetric transient simulation 


GASP 


Gulick et a!., 2006 
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Plate shock formation and collapse 
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Ground pressure 
spatial profiles 


Plate shock 
collapse 


— Quasi-Steady State (28 ms) 

Pc Ramp Down (54 ms) 

Plate Shock Collapse (61 ms) 

Complete Pc Ramp Down (88ms) 


Pc -1200 kPa 
0 deg cant 
h/de = 25 

- N 2 test gas 
-Axisymmetric 

- Transient 
-10 Hz pulsing 


FLUENT 


Plemmons, Mehta et al, 2008 
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Plate shock formation and collapse FLUENT 


Plemmons, Mehta et al, 2008 



MARS ATMOSPHERE ~ 700 Pa 


Pc - 1200 kPa 
0 deg cant 

Modeled as a 60 degree wedge to h/dg = 25 3-D transient simulations of full-scale 

reduce computational resources 10 Hz pulsing Phoenix piumes interacting at surface 

Equivaient piume gas 



3-D numericai simuiations show that ground pressure ioads are asymmetric and deveiop overpressures 
during rapid engine startup and shutdown 


Gulick et at, 2006 
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P/Pc-max {solid lines) 


Good agreement between experimental results and numerical simulations (CFD) 


Comparing spatial and temporal ground pressure profiles 

SPATIAL 


MARS ATMOSPHERE ~ 700 Pa 


TEMPORAL 


0.97 





CFD: 
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Dashed lines - surface pressure 
Solid lines - thruster inlet stagnation pressure 
(chamber pressure) 
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Red line - CFD 
Dots - Test Data 


Comparing plume 
shock structure 



(a) (b) (c) 

Shadowgraph FLUENT 


P,, - 1200 kPa 
0 deg cant 
h/de = 8.4 
10 Hz pulsing 

N 2 test gas 


GASP 


Gulick et at, 2006 
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Descent stage 

Baffle ^ 



t 

Baffle 


MLE 1 
2010MLE-1 


MLE 4 
2010MLE-4 


MLE6 ^ 
2010MLE-S 


MLE 7 

2010MLE-7 


MLE6 

2010MLE-6 


MLE 5 
2010MLE-5 


MLE2 

2010MLE-2 



Testbed 


Thruster 


Impingement 
Plate 




Variable 
Gfoiind $lop« 
Plats/Sand Bex 


- Va, scale MSL MLE nozzle 
-N2 test gas 
-steady operation 
-Mars atmosphere 
-NASA Ames Research Center 
Planetary Aeolian Laboratory 


gaiffle 


tmpji>gem 
Plate 
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- NASA OVERFLOW 2.1 


- 3-D time-marching implicit code 

- structured overset grid 

- Navier-Stokes eqns solved overfull domain and internal nozzle 

- SST turbulence model 

- compressibility correction 

- steady-state 

- frozen flow used for modeling rocket plume gases 
- 12 million cells 



0 0.2 0.4 0.6 0.8 1 1.2 ® Oz n.4 0.6 o.b i i.i 




Pc -1700 kPa 
22.5 deg cant 
h/de = 35 
N 2 test gas 
steady operation 


Repetitive 
overpressures not 
observed 
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Pc -1700 kPa 
22.5 deg cant 
h/de = 35 
N 2 test gas 
steady operation 


P IP 



Shows that the plume is 
collimated and leads to 
large pressure gradients at 
Mars atmospheric pressure 
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MARS ATMOSPHERE ~ 700 Pa 




P^-1.7e3 kP^ 
22.5 deg cant 
h/de --35 
Equivalent 
Plume gas 
steady-state 


fu sea e 
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NOTE; NoCIo mh. 



Plume is collimated 
and does not 
dissipate at large 
axial distances - in 
agreement with test 
data 
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OVERFLOW 


Pc ~ 1.7e3 kPa 
22.5 deg cant 
h/dc ~35 
N 2 test gas 
steady-state 
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MARS ATMOSPHERE ~ 700 Pa 
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OVERFLOW 
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CFD — Mach Contours 


Mars Ambient 
Pressure 


Earth Ambient 
Pressure 
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e = 4.50 


e- 0.02 


Mars - moderately underexpanded 
plumes lead to max ground 
pressure loads due to collimated 
plume structure and development 
of a small areal plate shock 

Earth - highly overexpanded 
plumes dissipate/no plate shock 
formation 

Moon - highly underexpanded 
plumes leads to a large areal plate 
shock - decreases ground 
pressure 
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Spatial Asymmetry 

SURFACE 
SHEAR STRESS 


MARS ATMOSPHERE ~ 700 Pa 



SURFACE PRESSURE 
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MACH CONTOUR 


P^~1200 kPa 
0 deg cant 
h/de = 25 
Steady-state 

N 2 test gas 

FLUENT 
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♦ Moderately underexpanded jets demonstrate: 

• collimated shock structures 

• large supersonic core lengths 

• plate shock dynamics 

• max pressure loads 


MISSION SUGGESTION: Due to highly 
complex plume impingement physics, 
accurate landing environments are 
needed for future planetary manned 
and robotic spaceflight missions 


♦ Plate shock dynamics leads to: 

• large pressure gradients 

• asymmetry 


j 


• overpressure 

♦ Ground pressure loads are highly sensitive to 

• jet expansion ratio 

• strouhal number 

• spacecraft altitude 

♦ Scaling laws show that cold plume gases can simulate ground pressure loads and interaction physics 
due to rocket plumes provided dynamic similarity is satisfied 

♦ How does this effect spacecraft landing? 

• Transient ground pressure loads translate to load perturbations at the spacecraft base which may lead to 
destabilizing moments (observed to a minor degree on the Phoenix spacecraft, Gulick et al, 2006) 

• Pressure loads can lead to extensive cratering and dust lifting which can destabilize the spacecraft upon 
touching down on the surface (observed at the Phoenix Landing Site, Mehta et al, 2011) 

• Dust lifting can erode important spacecraft sensors and science instrumentation (a concern for the MSL 
mission, Mehta et al, 2011) 

• Propulsion systems of small scale landers show maximum ground pressure loads at Mars atmosphere at 
relatively high altitudes (h/dg ~35) (a concern for the MSL mission) 

♦ Provided JPL with landing environments which they incorporated into their risk analysis models 
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-200 Hz (Inertial Measureiment Unit) IMU data and 10 Hz Radar data 



Current 

Research 

Investigation 
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Plumes interacted with the surface for less 
than 2 seconds (even less than predicted 
by landing simulations) 

Lift loss occurred at ~4.5 m and ground 
effect started around ~3.5 m 

Noticed a second bounce - could be the 
result of plume-surface interactions after 
Initial contact. 
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